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Einstein’s general theory of relativity predicts that an initially plane wave-front will curve because 
of gravity. This effect can now be measured using Very Long Baseline Interferometry (VLBI). A 
wave-front from a distant point source will curve as it passes the gravitational field of the Sun. We 
describe an idealised experiment to directly measure this curvature, using four VLBI stations on 
earth, separated by intercontinental distances. Expressed as a time delay, the size of the effect is a 
few hundred picoseconds and may be measureable with present technology. 
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It is now possible to do interferometry with intercon¬ 
tinental baselines (VLBI) and clocks of picosecond ac¬ 
curacy [1,2]. This technological advance of the last few 
years suggests a new test of general relativity. The test 
consists of directly measuring the curvature of a wave- 
front coming from a distant radio source. In the absence 
of general relativistic (GR) effects, the wave-front from 
the distant source would appear plane. General relativity 
predicts that the wave-front will curve because of gravi¬ 
tational effects. By means of four VLBI stations located 
on the earth, one can directly measure this curvature. 
The effect, whose size we estimate below, may be mea¬ 
surable with currently available VLBI techniques. We de¬ 
scribe an idealised version of the experiment in the hope 
of motivating VLBI astronomers to design more realistic 
experiments to detect the curvature of wave fronts. 

That gravity curves wave-fronts is already established 
from astronomical observations of gravitational lensing. 
Multiple images, caustics and “luminous arcs” prove that 
initially plane wave-fronts curve under the influence of 
gravity. However, all these gravitational lenses are out¬ 
side the solar system and we can only guess at and imper¬ 
fectly model their structure. Tests of general relativity 
within the solar system are far more under our control, 
since we are on familiar ground. We can model the lenses 
in detail, work out the predictions of the theory and quan¬ 
titatively confront theory with experiment. 

Two general relativistic effects in the propagation of 
light have already been accurately measured-the Shapiro 
time delay [3] and the bending of light [4,5]. It is worth 
noting that analogous effects exist even in special relativ¬ 
ity. The bending of light manifests itself as a change in 
the apparent direction of the source. In special relativity, 
uniform motion of the observer can result in aberration- 
an apparent change in the position of the source. Time 
delays too can be induced by uniform motion of the ob¬ 
server as in the Doppler effect. However curving of plane 
wave-fronts is a purely general relativistic effect that has 
no special relativistic analog. A plane wave in one in¬ 
ertial frame appears plane in all inertial frames. The 
experiment proposed in this paper consists of directly 
measuring the curvature of a wave-front, a purely GR 


effect with no special relativistic analog. 

This letter is organised as follows: we first describe 
the theory behind the proposed experimental test. We 
then describe the experiment in an idealised form and 
show how to extract the general relativistic effect from 
measured quantities. 

Consider a plane wave incident on an isolated static 
spherically symmetric body (the sun) of mass m and ra¬ 
dius R. We note that that e := 2 Gm/(c 2 R) = 2.10 -6 << 
1 and neglect terms of order e 2 . The effect of gravity can 
be calculated using the Schwarzschild metric. Neglecting 
terms of second and higher order in e, the Schwarzschild 
metric is given in standard co-ordinates by (we hence¬ 
forth set c = 1). 

ds 2 = (1 — p/r)dt 2 — (1 + p/r)dr 2 — r 2 dfl 2 , (1) 

where p = 2 Gm is the Schwarzschild radius of the body. 
We can write the metric g llv = + h^ v + 0(e 2 ), where 

„ is the General Relativistic perturbation of the metric 
and our notation 0(e 2 ) means that we neglect quantities 
which are second order and higher in e. The propagation 
of electromagnetic waves [6] can be described in the ge¬ 
ometrical optics approximation by the eikonal equation 

[7] 

= 0 

and the solution ip(x) is the phase of the wave, a scalar 
function of the general co-ordinate x^,p = 0,1, 2,3. Let 
us choose Minkowskian coordinates = (f, r) adapted 
to the flat metric Let lj be the frequency at infinity 
of the incident plane wave and let k^f be wave vector at 
infinity: kff = (w, 0, 0, — w), where we have chosen the 
direction of the incident plane wave along the positive z 
axis. In Minkowski space the function ip M (x) = kffx 11 
with rf 1 ' kf? = 0 solves the eikonal equation. To first 
Post-Minkowskian order, the solution is ip = ip M + cj), 
where (f> satisfies the differential equation 

p^kfd v <j> = 1 /2h^k™k ™. (2) 

Then (2) leads to 
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where f 
yields: 


u 2 + x 2 + y 2 and —> 

<t> = pui ln[(r — z) + z/2r\ 
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where r = \Jx 1 + y 2 + z 2 . For lines of sight close to the 
sun (x,y << z) it is more illuminating to use the simple 
approximate form 

(j)= {l±2)P^ [Hx 2 +y 2 )/{2z) \ (5) 

which we arrive at using a binomial expansion and drop¬ 
ping a constant independent of ( x,y,z ). In (5), we have 
introduced the PPN (parametrised post-Newtonian [5]) 
parameter 7 , which is equal to 1 in general relativity. 

The total phase of the wave is given by 

ip = + cp, (6) 

where <p is the small general relativistic correction, the 
Shapiro delay. The arrival of a wave-front at (x, y, z) is 
delayed by a time <p(x, y, z)/u> due to gravitational effects. 
Thus the gravitational field acts as a retarder, just as a 
glass slab does with light. There is a close mathematical 
relation between the Shapiro delay [3], the bending of 
light [4] and the curvature of wave-fronts which is the 
subject of this paper. Let us expand the eikonal ip in a 
Taylor expansion about an event x 3 : 

ip{x») = ip(x, 0 ) + (a^ - x%)d,j,ip\ Xo + 

1/2(2^ - x^)(x v - Xo)dnd v ip\ Xo + ... (7) 

The Shapiro delay measures ip(xo), the eikonal. The 
bending of light is a measurement of its first derivative 
kfj, = d^iplxo and the curvature effect described here mea¬ 
sures its second derivative. While these effects are re¬ 
lated mathematically, they are distinct physical effects 
and should therefore be separately measured and checked 
against theoretical predictions. In these three effects the 
gravitational held acts respectively as a retarder, a prism 
and a lens, which act by delaying, tilting and curving a 
wave-front. 

In the experiment to be described below, all observa¬ 
tion points are on the earth (see figure 1). As a result 
the numerical value of h the perturbation of the metric 
tensor is around 10 -8 . We will therefore neglect h /lv at 
the earth’s location and use the Minkowskian metric 77 
to raise and lower indices. Henceforth all dot products 
are formed using the hat Minkowskian metric. 

A wave-front is a three dimensional surface in four di¬ 
mensional space-time. For easy visualisation let us pick 
a slice t = 0. The wave-front is now a two dimensional 
surface described by 


—k.x + (p(x) = constant. 

In order to detect the curvature [8] of this wave-front 
one needs to sample at least four points on it. For, given 
three points (or less), one can always find a plane passing 
through all of them. The deviation of the fourth point 
from the plane passing through the first three gives a 
measure of the curvature of the wave-front. More sym¬ 
metrically, one can compute the volume of the tetrahe¬ 
dron with these four points as vertices. A non-zero vol¬ 
ume would imply that these four points do not he in a 
plane. This simple three dimensional argument provides 
the intuition we will use in describing the experiment. 
The principle of the experiment is that given the arrival 
events at three VLBI stations, one can model the inci¬ 
dent wave as a plane wave and absorb any deviation in a 
redefinition of the apparent direction of the source. How¬ 
ever, with the fourth station this freedom does not exist 
and one can measure genuine curvature effects [9]. Our 
method does not require an absolute determination of 
the direction of the source. 

Consider four VLBI stations l a , a = 0,1,2, 3 with in¬ 
tercontinental separations viewing a radio source along a 
line of sight passing near the sun. The VLBI experiment 
involves simultaneously observing a structureless point 
radio source from four telescopes. The signal voltage 
received at each antenna is beaten with an ultra-stable 
local oscillator, low pass filtered, digitized and recorded. 
The recorded data from two antennae are cross-correlated 
with a delay to find the ‘fringe’- the delay at which the 
cross-correlation peaks. Thus one can accurately locate 
the events at which a given wave-front arrives at the four 
antennae. 

Such VLBI techniques accurately measure the arrival 
times t a of a wave-front at each of these stations located 
at x a . Here x% = (t a . x a ) are the four four-vectors used to 
locate the four events, the arrival of a wave-front at each 
of the four stations. We will now describe how the mea¬ 
sured quantities x%, a = 0,1,2, 3 can be used to detect the 
curvature of the wave-front. Let us define = x% — Xq, 
the three baselines [10] connecting the reference station 
Xq to the other three stations. We write the components 
of Xjp as (X a ,Y a ,Z a ,T a ). Note that the baselines have 
spatial as well as temporal components. We now con¬ 
struct the following determinant from these twelve num¬ 
bers: 


S = det 


V! V 2 X 3 

Yi Y 2 Y 3 

(Z 2 -T 2 ) ( Z 3 -T 3 ) 


(9) 


If the wave-front were plane with wave-vector kjf = 
(o>,0, 0, — u), the last row in the determinant S would 
vanish and so would S. If the incident wave is plane, and 
has wave-vector = (ko,ki,k 2 ,k 3 ), where k\,k 2 are of 
order e, we find again that S vanishes: From k^k* 1 = 


2 



kl -k\-kl- kl = 0, we conclude that to first order in 
e, jfeo = k 3 . It then follows from 

k.X a = k 0 (T a - Z a ) - hXa - k 2 Y a = 0 (10) 

that the last row of (9) is a linear combination of the first 
two and so S vanishes. A non-vanishing S implies that 
the wave-front is not plane and therefore is a diagnostic 
for the curvature of the wave-front. 

Theoretically the measured arrival time difference T a 
between stations a and 0 can be written (to first order in 
e) as 

T a = Z a ~ T a (11) 

where Z a is the expected arrival time difference in the ab¬ 
sence of gravitational effects and r a = {^{Xa) — (f>{ xq))/lu. 
Plugging this into the expression for S we find that 

S = (nai + r 2 a 2 + r 3 a 3 ) (12) 

where aq = X 2 Y 3 — X 3 Y 2 (and cyclic) are the areas of the 
parallelograms formed by the spatial baselines projected 
on the x ~ y plane. 

Let us express r a using a Taylor expansion for <j>{x) 
keeping terms up to second order in x — Xo, as in (7). It 
is convenient to choose co-ordinates so that the reference 
event Xo is in the x — z plane, Xq = 0. Setting Xq = b, 
the impact parameter, we find that 

r a = (1 + 7 )p[XJb + (V a 2 - X 2 )/b 2 ] + ... (13) 

As expected, the first term, being linear in X a drops out 
of the expression for S and S is given by 

— t - Xl ) + a 2 {Y 2 - A 2 ) + a 3 {Y 2 - Af)] 

(14) 

A non-vanishing S signals curvature of the wave-front. 
We can express the effect in time units. Given all the 
measured quantities x% except t 3 , if the wave-front were 
plane we would expect that the arrival time of the wave- 
front at station three would be t 3 where t 3 satisfies the 
equation 

(Zi — Ti)ai + ( Z 2 — T 2 )a 2 + (Z 3 — t 3 + t 3 )a 3 = 0, (15) 

This expectation would be false, because S does not van¬ 
ish and we have 

( Z\ — Ti)ai + (Z 2 — T 2 )a 2 + (Z 3 — t 3 + to)a 3 = S (16) 

Subtracting (15) from (16) we find that t 3 would differ 
from t 3 by an amount A t 3 = t 3 — t 3 given by 

At 3 = S/a 3 (17) 

The difference A< 3 between the measured arrival time t 3 
and the naive expectation t 3 would detect the curvature 


of the wave-front. The size of the effect is easily estimated 
from (12,17). The main purpose of this letter is to draw 
attention to the fact that this effect may be measureable 
with present VLBI technology. Expressed in time units, 
it is of order 

10~ 5 |A| 2 /6 2 s, (18) 

where |A| is the typical size of the baseline (the intercon¬ 
tinental separation of the telescopes) and b the impact 
parameter. If one chooses the impact parameter to be 
three solar radii, (b = 3 Rq), the effect is 100 ps, which 
may be measureable. At grazing incidence, ib = Rq), the 
effect is as large as a nanosecond. However, lines of sight 
close to the Sun suffer from the problem of noise due to 
the Solar ionosphere. As one chooses larger impact pa¬ 
rameters, this noise is reduced, but so also is the effect 
of interest (which decreases as the square of the impact 
parameter b). The optimal choice of impact parameter 
may be best determined by trial and error. 

The ideal source for this experiment would be a strong, 
distant point source. A number of such sources have been 
already identified in the ICRF (International Celestial 
Reference Frame) catalog [11]. These sources are densely 
distributed over the sky with an average separation of a 
few degrees and can be detected using integration times 
of just a few minutes. Apart from the non-dispersive 
gravitational effects of interest there are some dispersive 
non-gravitational effects such as due to ionospheric fluc¬ 
tuations. Such phases can be removed by dual or multi¬ 
frequency observations, which is a standard technique in 
VLBI. There are also non-dispersive, non-gravitational 
effects such as due to the troposphere. These can be re¬ 
moved by the technique of phase referencing [1] between 
sources, provided the sources are within a few degrees of 
each other. 

As we mentioned earlier, the effect described here is 
contained in the Shapiro effect, just as the bending of 
light is. Since the Shapiro time delays are routinely ac¬ 
counted for in VLBI observations, the proposed effect 
has already been implicitly measured. The purpose of 
this letter is to motivate observers to explicitly measure 
this effect that general relativity predicts. Although the 
Shapiro effect has been measured with high precision un¬ 
der varied conditions [12], this does not constitute a di¬ 
rect measurement of the curvature of a wave front, since 
one does not sample four points on the same wavefront. 
We suggest here that a the curvature of the wave front 
is worth measuring directly. Indeed, the mean curvature 
of the wave-front is related via Raychaudhuri’s equation 
[13] to the integral of the Ricci curvature along the line 
of sight. In Einstein’s theory, we expect the wavefront to 
be a minimal surface (zero mean curvature). This pre¬ 
diction could be explicitly checked. To conclude, we have 
proposed an idealised experiment to directly measure the 
curvature of the wave-front from a distant radio source. 
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We hope to motivate VLBI astronomers to design a more 
realistic version of the experiment to measure the rela¬ 
tivistic curvature of wave fronts. 
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FIG. 1. Figure 1 shows disposition of the radio telescopes, 
all on the same hemisphere as the source SI, which is in com¬ 
mon view of all the telescopes. The source S2 about three 
degrees away from SI is used for phase referencing. The fig¬ 
ure is not to scale. 
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